Crossover to Potential Energy Landscape Dominated Dynamics 
in a Model Glass-forming Liquid 
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An equilibrated model glass-forming liquid is studied by mapping successive configurations produced 
by molecular dynamics simulation onto a time series of inherent structures (local minima in the 
potential energy). Using this "inherent dynamics" approach we find direct numerical evidence for 
the long held view that below a crossover temperature, T x , the liquid's dynamics can be separated 
into (i) vibrations around inherent structures and (ii) transitions between inherent structures (M. 
Goldstein, J. Chem. Phys. 51, 3728 (1969)), i.e., the dynamics become "dominated" by the potential 
energy landscape. In agreement with previous proposals, we find that T x is within the vicinity of the 
mode-coupling critical temperature T c . We further find that at the lowest temperature simulated 
(close to T x ), transitions between inherent structures involve cooperative, string like rearrangements 
of groups of particles moving distances substantially smaller than the average interparticle distance. 

PACS numbers: 61.20.Lc, 61.20.Ja, 63.50.+X, 64.70.Pf 



I. INTRODUCTION 

Dynamical behavior of many physical and biological 
systems JlljJ] can be considered in terms of the transient 
localization of the system in basins of potential energy, 
and transitions between basins. In particular, this ap- 
proach has received much attention in studies of slow 
dynamics and the glass transition in supercooled liquids. 
Here, the strong temperature dependence of transport 
properties such as the diffusion coefficient and viscosity, 
and the possible existence of a thermodynamic transi- 
tion underlying the laboratory glass transition, have been 
sought to be understood in terms of the properties of 
the liquid's potential energy (or free energy) surface, or 
"landscape" as it is commonly called JHp|jF^lf|| . 

For a system composed of N atoms, the potential en- 
ergy surface is simply the system's potential energy plot- 
ted as a function of the 3N particle coordinates in a 
3 N + 1 dimensional space || . The potential energy sur- 
face contains a large number of local minima, termed 
"inherent structures" by Stillinger and Weber [[|. Each 
inherent structure is surrounded by a "basin" , which is 
defined such that a local minimization of the potential en- 
ergy maps any point in the basin to the inherent structure 
contained within it. The time evolution of a liquid may 
be viewed as the motion of a point on the potential energy 
surface, and thus as a succession of transitions from one 
basin to another. These transitions are expected to occur 
differently as the temperature T is varied. In particular, 
Goldstein argued j|] that below a crossover temperature, 
T x , where the shear relaxation time is ~ 10~ 9 seconds, 
relaxation is governed by thermally activated crossings 
of potential energy barriers. The presence of significant 
energy barriers below T x suggests a clear separation of 



short-time (vibrational) relaxation within potential en- 
ergy basins from long-time relaxation due to transitions 
between basins. 

A complementary approach to the dynamics of super- 
cooled liquids is provided by the mode coupling theory 
(MCT) 0. The simplest (so-called "ideal") version of 
this theory predicts a power-law divergence of relaxation 
times and the inverse diffusion coefficient, at a critical 
temperature T c . Although a power law provides a reason- 
able description of the temperature dependence of these 
quantities above T c in both real and simulated systems, 
power law behavior breaks down for T ks T c , i.e. the 
predicted singularity at T c is not observed. This devia- 
tion is attributed to the presence of "hopping" motion 
as a mechanism of relaxation, which is not included in 
ideal MCT Jl7[ ]. Consequently, T c is usually estimated 
by fitting a power law to a relaxation time, taking into 
account that this fit is expected to break down close to 
(and below) T c . 

It was noted by Angell that experimentally it is of- 
ten found that the shear relaxation time is on the order 
of 10~ 9 seconds at the estimated T c , leading to the argu- 
ment that T x pa T c (See also Ref. ]18|] ). The presence of a 
low temperature regime where barrier crossings dominate 
the dynamics, and the correspondence of the crossover to 
that regime with the mode coupling critical temperature 
T c , has also been discussed in the context of mean field 
theories of certain spin glass models 

The existence of a crossover temperature and corre- 
sponding separation of the dynamics can be directly 
tested with computer simulations, using the concept of 
inherent structures. In this paper, we map the dynamical 
evolution of an equilibrated model liquid to a time series 
of inherent structures for a range of temperatures. In this 
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way, we test the extent to which short-time "intra-basin" 
relaxation is separable from long-time "inter-basin" re- 
laxation. Our results demonstrate that this separation 
becomes valid as the system is cooled, and we estimate 
the crossover temperature T x to be close to the estimated 
value of T c . 

II. INHERENT DYNAMICS 



In this paper, we quantitatively compare F s (q,t) 
and Fg(q,t) to test whether the dynamics of a binary 
Lennard-Jones mixture can be separated into vibrations 
around, and transitions between inherent structures. If 
so, then F^(q,t) describes the relaxation of the liquid as 
described by F a (q,t), but with the effect of the vibrations 
removed. We show that this scenario becomes true below 
a crossover temperature, T x , which is close to the lowest 
temperature simulated in the present work. 



In this section we describe the details of our approach, 
which is sketched in Fig. [I]. After equilibration at a 
given thermodynamic state point, a discrete time series 
of configurations, R(i), is produced by standard molec- 
ular dynamics (MD) simulation. Each of the configura- 
tions R(i) is then mapped to its corresponding inherent 
structure, R/(i), by locally minimizing the potential en- 
ergy in configuration space. We refer to this procedure as 
a "quench". After quenching the configurations in R(£), 
we have two "parallel" time series of configurations, R(i) 
and R / (i). The time series R(£) defines the "true dynam- 
ics" , which is simply the usual (Newtonian) MD dynam- 
ics. In an analogous way, the time series R J (t) defines 
the "inherent dynamics" . If a function quantifying some 
aspect of the true dynamics is denoted by f(R(t)), then 
the corresponding function, /(R I (i)), of the inherent dy- 
namics is calculated in exactly the same way, except using 
the time series of inherent structures. For example, the 
self intermediate scattering function, F s (q,t), and the in- 
herent self intermediate scattering function, F*(q,i), are 
defined by 



F s (q,t) = (cosq-(r,(i)-r,(0))} , 
F/( 9 ,i) = <cosq-(rj(t)-rj(0))) 



and 



(2.1) 
(2.2) 



where r£(t) is the position of the jth particle in the in- 
herent structure R J (t) and (...) denotes an average over 
j and the time origin. 
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FIG. 1. Schematic describing the principle of the "inherent 
dynamics" approach. Successive configurations of the equili- 
brated liquid R(t) are quenched to produce their correspond- 
ing inherent structures R 7 (t). Successive inherent structures 
form a time series which we use to calculate the inherent self 
intermediate scattering function F%(q, t) (Eq. |2.2| ). More gen- 
erally, the inherent counterpart of any equilibrium quantity 
may be calculated in this fashion. 



III. RESULTS 

In the following we present results from molecular dy- 
namics simulations of a binary Lennard-Jones mixture 
in three dimensions, equilibrated at eight different tem- 
peratures. The model used for the present simulations 
is described in Ref. |p2fl . The system contains 251 par- 
ticles of type A and 249 particles of type B interacting 
via a binary Lennard-Jones potential with parameters 
<Jbb/<?aa = 5/6, (T AB = (f aa + <Jbb)/2, and e AA = 
cab = £ bb ■ The masses are given by mj/mj = 1/2. 
The length of the sample is L — 7.28<taa and the po- 
tential was cut and shifted at 2.ba a p. All quantities are 
reported in reduced units: T in units of caa, lengths in 
units of <jaa and time in units of r = (mscr^ yi /48e) 1 ' 2 
(this was misprinted in |22|). Adopting "Argon units" 



leads to a A A = 3.4A, e/k B = 120K, and t = 3 x l(T 13 s. 
The simulations were performed in the NVE ensemble 
using the leap-frog algorithm with a timestep of O.Olr, 
at constant reduced density, p = 1.296. The quenching 
was performed using the conjugate gradient method p3[ . 

We first briefly describe aspects of the true dynamics 
that demonstrate a qualitative change occuring in the 
temperature range investigated. 

In Fig. H we show the quantity Airr 2 G s A{r, ti), which 
is the distribution of displacements 0] of particles of 
type A during the time interval t\. We define t\ as 
the time where the mean square displacement is unity, 
(r 2 (ti))A = 1. At all temperatures the dynamics be- 
come diffusive {{r 2 {t))A °c t) for t > t\ (see inset), i.e., 
t\ marks the onset of diffusivity. At the highest tem- 
peratures, 4:irr 2 G s A(r, fi) agrees well with the Gaussian 
approximation [thick curve, G s A{r,ti) cx exp(— 3r 2 /2)]. 
As T is lowered, the distribution of particle displace- 
ments deviates from the Gaussian approximation, and 
a shoulder develops at the average interparticle distance 
(r « 1.0 in the adopted units), which at T=0.59 becomes 
a well-defined second peak. The second peak, observed 
also in other model liquids at low temperatures, indicates 
p5| , p6]| single particle "hopping" (see Fig. ||a): particles 
stay relatively localized for a period of time (first peak), 
and then move approximately one interparticle distance, 
where they again become localized (second peak). Thus 
we see from Fig. |^ that as we approach our lowest simu- 
lated temperature T = 0.59, there is a qualitative change 
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from dynamics well described by a Gaussian distribution 
to dynamics dominated by hopping processes. 
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FIG. 2. Distribution of particle displacements for the 
A (large) particles, 4ivr 2 GsA(r , ti), where ti is denned by 
(r 2 (ti))A = 1 (see inset). At high T the Gaussian approx- 
imation (thick curve) is reasonable, whereas at the lowest T 
a second peak is present, indicating single particle hopping. 
Inset: Mean square displacement of the A particles, (r {t))A- 
Similar behavior is found for the B (small) particles. 
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FIG. 3. (a) Trajectory of a particle at T = 0.59. The 
elapsed time is At = 160r (the typical "vibration" time is 
ss It). At this temperature the dynamics is dominated by 
"hopping": particles stay relatively localized for many time 
steps and then move approximately one interparticle distance, 
where they again become localized, (b) Applying the inherent 
dynamics approach to the trajectory above. The 1600 con- 
figurations used to generate the (true) trajectory in (a) were 
quenched, and the positions of the particle in the resulting in- 
herent structures are here plotted and connected by straight 
fines. 

In Fig. [|b the inherent dynamics approach is applied 
to the true trajectory seen in Fig. ^a. The resulting "in- 
herent trajectory" consists of the positions of the particle 



in 1600 successive quenched configurations. The quench- 
ing procedure is seen to remove the vibrational motion 
from the true trajectory. The inherent trajectory will be 
discussed in more detail in section IV. 




FIG. 4. (a) F s a (q, t) plotted versus t on log-scale for 
q — 7.5 at the same temperatures as in figure ^[ Data 
points are connected by straight lines. Dashed lines are fits to 
fit) = f c exp( — (t/Ta) 13 ). (b) Fs A (q, t), otherwise as above. In 
both (a) and (b), the fitting was performed for t > 10 for the 
two highest temperatures and for t > 30 otherwise. Similar 
behavior is found for the B particles. 

We now compare the true self intermediate scattering 
function, F s (q,t), with its inherent counterpart Fj(q,t). 
Fig. |I|a shows the self intermediate scattering function for 
the A particles, F s A{q, t), at q = 7.5 corresponding to the 
position of the primary peak in the static structure factor 
for the A- A correlation. For each temperature F s (q,t) 
was calculated from approximately 2000 configurations 
(depending on temperature). As T decreases, F s A(q,t) 
is found to display the typical two-step relaxation, where 
the short time decay is attributed to vibrational relax- 
ation (or "dephasing", see Ref. of particles within 
cages formed by neighboring particles ]27|-p9| . The long 
time, or a-relaxation is separated from the short time 
regime by a plateau indicating transient localization, or 
"caging" of particles, and is generally observed to follow 
a stretched exponential form. 

The self part of the inherent intermediate scattering 
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function for the A particles, F* A (q,t) at q=7.5, is shown 
m Fig. |b. This was calculated by quenching each config- 
uration used in Fig. pi, and then applying the same data 
analysis program on the resulting time series of inherent 
structures. As expected, the plateau disappears in the in- 
herent dynamics, as previously shown also for the inher- 
ent mean-square displacement |E| . At all T we find that 
the long-time behavior of both F s A(q,t) and F^ A (q,t) is 
well described by stretched exponentials (dashed lines). 
As a result, we can quantitatively compare the long time 
relaxation of F s A(q,t) and F^ A (q,t), by comparing the 
fitting parameters {r Q ,/3, / c } of the stretched exponen- 
tials f(t) = f c exp(-(t/T a f). 

If the true dynamics can be separated into vibrations 
around and transitions between inherent structures, how 
do we expect the fitting parameters for the inherent self 
intermediate scattering function, {r^, (3 1 , f£ } to be re- 
lated to the fitting parameters for the true self inter- 
mediate scattering function, {r Q ,/3, / c }? To answer this 
question, we assume that the initial relaxation in F s (q, t) 
is due to vibrations (as widely accepted [ plp^ -[29|l). If 
this is the case, then we expect the quenching procedure 
to remove the initial relaxation (since it removes the vi- 
brations), which means that F^(q,t) can be thought of 
as F s (q,t) with the initial relaxation removed^. This in 
turn means that -F 1 / (q, t) should be identical to the long 
time relaxation of F s (q,t), but rescaled to start at unity: 
{ri,/3 / ,/ c / } = {r Q ,/3,l}. 

The fitting parameters used for fitting stretched ex- 
ponentials to F sA (q,t) (Fig. |a) and F' A (q,t) (Fig. |b) 
are shown in Fig. ^ (a) relaxation times, r Q and t^, (b) 
stretching parameters, f3 and and (c) non-ergodicity 
parameters, f c and f 1 . We also show in Fig. |5|a the 
fit of the asymptotic mode coupling prediction r Q oc 
(T - T c )-t, from which we find T c = 0.592 ± 0.006 and 
7 = 1.41 ± 0.07. The fitting was done without the low- 
est temperature, where hopping is clearly present in the 
system (see Fig. |]), since this type of particle motion 
is not included in the ideal mode coupling theory. Ex- 
cluding the two lowest T gives a fit which is consistent 
with the one presented here; including all temperatures 
gives a considerably worse fit. Applying the same pro- 
cedure to the inverse diffusion coefficient, Z? _1 (T), gives 
T c = 0.574 ± 0.005 and 7 = 1.40 ± 0.09 (data not shown). 



*If vibration can be separated from transitions between 
inherent structures we may write for the x-displacement: 
Ax = Ax V ib + Axi n h where the two terms are statistically 
uncorrelated. Thus, [using an exponential instead of cosine 
in Eqs. (2.1) and (2.2)] we find that the self intermediate scat- 
tering function is a product of a term relating to vibrations 
and one relating to transitions between inherent structures. 
At long times the former becomes time-independent, converg- 
ing to the non-ergodicity parameter. 
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FIG. 5. Parameters describing the fit of F s A(q = 7.5, t) 
and F^ A {q = 7.5, t) to stretched exponentials from Fig. 3 (a) 
and (b), respectively, (a) Relaxation times r a and vs. 
T. The solid line is a fit to r a oc (T — T c ) ' excluding the 
lowest T in the fitting (see text). INSET: ri vs. r Q . (b) 
Stretching parameters f3 and /3 1 vs. T. INSET: /3 1 vs. (3. (c) 
Non-ergodicity parameters f c and fl vs. T. Error bars are es- 
timated from deviations between three independent samples. 
Similar behavior is found for the B particles. 

Also shown in Fig. |5| as insets are vs. t q and (3 1 
vs. (3. Within the error bars we find that r a and are 
identical at all temperatures. At the highest tempera- 
tures (3 is poorly defined since there is no well-defined 
plateau in F s A(q,t). Consequently it is difficult to com- 
pare /3 and (3 1 at high T, but we find that they become 
identical (within the error bars) at low T. Thus at low 
temperatures our results confirms the expectation that 
the inherent dynamics is simply a coarse-graining of the 
true dynamics, i.e., that {t^, ft 1 } = {r Q , /?}. On the other 
hand, the non-ergodicity parameters f c and (Fig. ||c) 
are strikingly different. While f c is roughly independent 
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of T, f£ increases towards unity as T approaches our low- 
est temperature. The fact that we observe a temperature 
dependence of f[ approaching unity as T approaches our 
lowest temperature T = 0.59, leads us to conclude that 
this is close to the crossover temperature, T x . We note 
that Goldstein's estimate of shear relaxation times at T x 
(1CT 9 seconds) in our LJ units corresponds to 3 x 10 3 , 
which is the same order of magnitude as r Q in the tem- 
perature range where f c approaches unity. 

Below T x the inherent dynamics can be thought of as 
the true dynamics with the effect of the vibration re- 
moved, as shown above. How should the inherent dy- 
namics be interpreted above T x l In Fig. ^b the short 
time relaxation of the inherent self intermediate scatter- 
ing function at high temperatures is seen to be approx- 
imately logarithmic in time. This is an artificial relax- 
ation introduced by applying the quenching procedure 
at a temperature where the dynamics is not separated 
into vibrations around, and transitions between inher- 
ent structures, i.e. the quenching procedure is doing 
more than simply removing the vibrations around inher- 
ent structures. Presumably the inherent dynamics above 
T x contains information about the underlying potential 
energy landscape. At the present, however, we do not 
know how to interpret this, and we do not have an expla- 
nation as to why the (artificial) initial relaxation appears 
to be logarithmic at high temperatures. 

We now proceed to discuss Angell's proposal, that 
T x « T c . We find that both estimated values for T c 
[0.592±0.005 from r a (T) and 0.574±0.005 from D~ l {T)} 
are in the temperature range where f£ is approaching 
unity. We note that in the system investigated here two 
of the asymptotic predictions of the ideal mode coupling 
theory do not hold; t q and D~ l have different temper- 
ature dependence and we do not find time-temperature 
super-position of the a-part of the self intermediate scat- 
tering function. However, the argument given by Angell 
jjj (and Sokolov |lq| ) only relates to T c as the temper- 
ature where power-law fits to experimental data tend to 
break down, i.e. the "usage" of MCT in this argument is 
similar to the way we have estimated T c in Fig. [|a, and 
does not require, e.g., time-temperature superposition. 



IV. TRANSITIONS BETWEEN INHERENT 
STRUCTURES 

As shown in the previous section, separation of the dy- 
namics into vibrations around and transitions between 
(the basin of attraction of) inherent structures becomes 
possible as T approaches T x , which is close to our low- 
est simulated temperature T=0.59. At this temperature, 
it therefore becomes meaningful to examine the details 
of the transitions between successive inherent structures. 
We identify such transitions by quenching the MD config- 
urations every O.lr (i.e. every 10 MD-steps) and looking 



for signatures of the system undergoing a transition from 
one inherent structure to another. We have considered 2 
such signatures: i) We monitor the inherent structure en- 
ergy E 1 ^) as a function of time, as shown in Fig. ||a. ii) 
We monitor the distance in configuration space AR 1 (t) 
between two successive quenched configurations Q (Fig. 
^b), where 



AR^t) = |R/(t + 0.1) -R/(i)| 



N 



£(r£(t + 0.1)-rj(i))' 



(4.1) 
(4.2) 



Each jump in E 1 (t) corresponds to a peak in AR 1 ^), 
indicating a transition to a new inherent structure. In 
the (rare) event where a transition occurs between two 
inherent structures with the same energy, AR 1 (t) will 
still exhibit a peak even in the absence of a jump in E 1 (t) , 
and for this reason we use AR 1 (t) to identify transitions. 
The condition AR 1 (t) > 0.1 was found to be a sufficient 
threshold for this purpose. When evidence of a transition 
was found in a time interval At = O.lr, this time interval 
was divided into 10 subintervals of At — O.Olr and the 
procedure described above was repeated. 
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FIG. 6. Identifying transitions between inherent struc- 
tures, (a) The inherent structure energy E 1 (t) vs. time, (b) 
AR 1 (Eq. ^l]) vs. time. A transition between (the basin of 
attraction of) two inherent structures is indicated by a jump 
in E 1 (t) and a corresponding peak in AR 1 (t). 

For each transition, we monitor the difference between 
the particle positions in the two successive inherent struc- 
tures. The distribution p(r) of all such particle "displace- 
ments" averaged over the 12000 transitions we have iden- 
tified is shown in Fig. ^. While many particles move only 
a small distance (r < 0.2) during a transition from one in- 
herent structure to the next, a number of particles move 
farther, and in particular, we find that the distribution 
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for r > 0.2 is to a good approximation exponential. The 
dotted curve is a fit to a power-law with exponent —5/2, 
which is a prediction from linear elasticity theory [ 61] , 
describing the displacements of particles in the surround- 
ings of a local rearrangement "event" . This power-law fit 
does not look very convincing by itself, but we note that 
the exponent was not treated as a fitting parameter (i.e. 
only the prefactor was fitted), and the power-law must 
break down for small displacements, since these corre- 
spond to distances far away from the local event, and are 
thus not present in our relatively small sample. From the 
change in behavior of p(r) at r w 0.2, it is reasonable to 
think of particles with displacements larger than 0.2 as 
those taking part in the local event, and the rest of the 
particles as merely "adjusting" to the local eventfj]. Using 
this definition it is found that on average approximately 
10 particles participate in an event. 
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FIG. 7. Distribution of particle displacements during 
transitions between successive inherent structures at T=0.59. 
The integral of the distribution is normalized to be the num- 
ber of particles N = 500. Full curve is a fit to an exponential, 
for 0.3 < r < 1.0. The dotted curve is a fit to a power-law 
with (fixed) exponent -5/2 [fl}, for 0.1 < r < 0.2 

Fig. ^ has two important consequences with regards 
to points discussed earlier in this paper. The first point 
relates to the single particle hopping indicated by the 
secondary peak in 47rr 2 G s (r, t) (Fig. |J) at low temper- 
atures. A common interpretation of the single particle 
hopping is that the jump of a particle from one "local- 
ized state" (first peak) to the second localized state (sec- 
ondary peak), corresponds to the transition of the sys- 
tem over an energy barrier from one inherent structure 
to the next. If such a transition typically occurs over 
a single energy barrier, i.e. without any new inherent 
structures between the two states, we would expect to 
find a preference for displacements of one average inter- 



^Note however, that our data does not imply what is cause 
and what is effect, or even if such a distinction is meaningful. 



particle distance (r w 1) in Fig. That this is not the 
case demonstrates that the hopping indicated by the sec- 
ondary peak in 4-7rr 2 G s (r, t) at low temperatures is not 
due to transitions over single energy barriers. Instead, as 
seen in the inherent trajectory in Fig. |^, the jump occurs 
via a number of "intermediate" inherent structures. 

The second important consequence of Fig. ^| is that 
particles in the surroundings of a local event are dis- 
placed by small distances. This kind of motion is diffi- 
cult to detect in the true dynamics, since it is dominated 
by the thermal vibrations. Presumably this kind of mo- 
tion is the reason why the inherent trajectory in Fig. || 
shows small displacements 0.2(7^), even when the 
corresponding true trajectory seems to oscillate around 
the same position: When a transition between inherent 
structures involving significant particle rearrangements 
in the surroundings occurs, the particle starts vibrating 
around a position that is slightly displaced, and a corre- 
sponding small displacement of the inherent trajectory is 
seen. This view of the dynamics is also consistent with 
the fact that the first peak in the inherent counterpart 
of 4irr 2 G s (r, t) (not shown, see Ref. p^] ) is not a delta 
function in r = 0. 




FIG. 8. Before (light) and after (dark) one typical tran- 
sition, all the particles which move a distance greater than 
0.2cr aa . Particles are shown with a diameter of <jaa- Note 
that most particles move considerably less than oaa (com- 
pare Fig. pj). The cooperative, string-like nature of the par- 
ticle motions during the inter-basin transition can be clearly 
seen. 

By observing, for a number of transitions, the positions 
of all particles that moved a distance greater than 0.2 
during a transition, we find these particles to be clustered 
together in "strings" , as shown in Fig. ||. Typically, one 
transition appears to involve just one string-like cluster. 
Detailed investigations of the transition events will be 
presented in a separate publication. Here we simply note 
that string-like particle motion has been observed also in 
the true dynamics above T c in a similar binary Lennard- 
Jones mixture []32"f . These strings are found on long time 
scales and involve particles moving approximately one 
inter-particle distance, and are thus different from, but 
presumably related to, the strings found in the present 
work. 
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V. CONCLUSIONS 

We have investigated the dynamics of a model glass- 
forming liquid in terms of its potential energy landscape 
by "quenching" a time series of MD configurations to 
a corresponding time series of inherent structures. In 
this way we have provided numerical evidence for the 
conjecture, originally made by Goldstein 30 years ago in 
this journal Q, that below a crossover temperature T x 
the dynamics of the liquid can be separated into vibra- 
tions around and transitions between inherent structures. 
Specifically, by comparing the self intermediate scatter- 
ing function F s (q, t) with its inherent counterpart F^(q, t) 
we presented evidence for the existence of T x . It is per- 
haps not surprising that the dynamics of a liquid becomes 
dominated by the structure of the potential energy land- 
scape at sufficiently low temperatures. What we have 
done here using the concept of inherent dynamics, is to 
provide direct numerical evidence for this, and we have 
shown that this regime can be reached by equilibrium 
molecular dynamics (for the particular system investi- 
gated here). To our knowledge this is the first time such 
evidence has been presented. 

In agreement with previous proposals j?],|l8],[l9) we find 
T x k. T c , where T c is estimated from a power-law fit 
to r a . This is also the temperature range where single 
particle hopping starts to dominate the dynamics, and 
r a becomes on the order of 10 -9 seconds (Goldstein's 
estimate of the shear relaxation time at T x ). 

The fact that we have been able to cool the system, 
under equilibrium conditions, to temperatures where the 
separation between vibrations around inherent struc- 
tures and transitions between these is (almost) complete, 
means that it becomes meaningful to study the individ- 
ual transitions over energy barriers, since the transitions 
in this regime dominate the dynamics. Our two key find- 
ings with regards to the individual transitions between 
inherent structures are i) single particle displacements 
during transitions show no preference for displacements 
on the order of the inter-particle distance, showing that 
the single particle hopping indicated in 4nr 2 G s (r, t) at 
low T (Fig. ||) does not correspond to transitions of the 
system over single energy barriers; and ii) particle dis- 
placements during transitions are spatially correlated (in 
"strings"). 
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